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Abstract 

The variational Monte Carlo method is used to find the ground state of six quarks confined 
to a cavity of diameter R c , interacting via an assumed non-relativistic constituent quark model 
(CQM) Hamiltonian. We use a flux-tube model augmented with one-gluon and one-pion exchange 
interactions, which has been successful in describing single hadron spectra. The variational wave 
function is written as a product of three-quark nucleon states with correlations between quarks in 
different nucleons. We study the role of quark exchange effects by allowing flux-tube configuration 
mixing. An accurate six-body variational wave function is obtained. It has only ~ 13% rms 
fluctuation in the total energy and yields a standard deviation of < .1%; small enough to be useful 
in discerning nuclear interaction effects from the large rest mass of the two nucleons. Results are 
presented for three values of the cavity diameter, R c = 2, 4, and 6 fm. They indicate that the 
flux-tube model Hamiltonian with gluon and pion exchange requires revisions in order to obtain 
agreement with the energies estimated from realistic two- nucleon interactions. We calculate the 
two-quark probability distribution functions and show how they may be used to study and adjust 
the model Hamiltonian. 
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I. INTRODUCTION 



The constituent quark model (CQM) has been useful to describe the spectroscopy [1], 0, |3| 
and the decay ||, |5|, [J of baryons and mesons. It retains only the constituent quark (CQ) 
degrees-of-freedom; the effects of all other degrees-of-freedom are subsumed into potentials 
dependent upon CQ positions, spins, flavors and colors. It assumes that a Hamiltonian 
describing interacting quarks can provide a useful description of the low-energy properties 
of hadrons. 

It is natural to attempt to extend the CQM to describe the two-baryon states made up 
of six CQ. The only such bound state known is the deuteron; all other known two-baryon 
states are unbound. The scattering data provides information on the interaction between 
baryons which is used to construct realistic models of the two-baryon potential [0, |]. § ■ 



Many authors [[TO], [11], |12| have attempted to calculate properties of two-baryon states, 
including the interaction between nucleons, using the CQM. In these studies some approxi- 
mation scheme is used to avoid the calculation of six-body eigenstates of the assumed CQM 
Hamiltonian. In the past few years there have been significant advances in the variational 
(VMC) and Green's function Monte Carlo (GFMC) methods to calculate the eigenstates of 



up to eight interacting nucleons [13, 14]. In these quantum Monte Carlo (QMC) methods 
a good approximate solution is first obtained by the VMC method from which the exact 
eigenstate is projected with the GFMC. 

In the present work we attempt to find the ground state of two-nucleons confined to a 
spherical cavity of diameter R c in the center-of-mass frame using the CQM. In the interior 
of the cavity the six CQ wave function is determined with VMC. Near the edge of the cavity 
we assume that it factorizes into a two-nucleon wave function required to be zero at the 
internucleon distance R c . 

This problem may be easily solved assuming that the ground state can be described as 
that of two interacting nucleons. The deuteron has total spin-isopin, S,T = 1,0, and in 
local, realistic models such as Argonne v\$ [|J, the nuclear interaction in the deuteron is 
given by a sum of central, tensor, and spin-orbit potentials: v^(r), t^(r) and vff(r). The 
two-nucleon Schrodinger equation: 

h 2 d 2 



2 -w(r) +v c d {r)u{r) + V8^(r)w(r) = E(R c )u{r) (1) 
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rriN \dr 2 r 2 



w{r) + ( v c (r) _ 2v t d{r) _ 3 ^ (r) ) ^ (r) 

+V8v t d {r)u{r) = E{R c )w{r) (2) 



where is the average mass of the neutron and proton, and u(r) and w(r), the 3 Si and 
! wave functions, can be easily solved with the boundary conditions u(R c ) = w(R c ) = 0. 
The results obtained with the v is model are shown in Table | for R c = 2, 4, and 6 fm. They 
presumably have small dependence on the model NN interaction because all models have 
mostly the one-pion exchange potential (OPEP), v^ N at r > 2 fm in T = states. Note 
that E{R C — > oo) is just the deuteron binding energy of —2.24 MeV, for the isoscalar part 
of Argonne V\%. 

The E(R C ) = (T N ) + (v^n), has a large cancellation between the nucleon kinetic energy 
(Tat) and the negative two-nucleon interaction energy (vnn), typical of nuclear systems [JT3 . 
In Table | we also list the energy, 

Em{Rc) = ^RV (3) 
of two non-interacting nucleons in a cavity of diameter R c and the total effect SE emp (R c ) = 
E(R C ) — E NI (R C ) of strong interactions, per nucleon. Our long-range aim is to calculate 
E(R C ), and thus SE(R C ), with computational errors of < 1 MeV using CQM, and compare 
with the SE emp (R c ), calculated from Argonne vi$ potential, which reproduces the known 
iViV data. These calculations can then be used to refine the CQM Hamiltonian and study 
the behavior of the wave function in the region where the six CQ are close together, i.e. 
the region in which the quark distributions of the two-nucleons overlap. Eventually such 
calculations may be useful to study A — N and £ — iV interactions for which the scattering 
data is limited. 

The first main concern in this approach is that the total energy, Ey, in the CQM includes 
the rest mass of the nucleons, and is therefore of order 2mjv ~ 2000 MeV. In order to calculate 
this energy with an accuracy of ~ 1 MeV the variance of the local energy, 

= <(R)J/»y(R) 
*;,(R)*v-(R) 

where R = {ri, T2, ■ ■ ■ , r$} is the six CQ configuration and \fV(R) is the variational wave 
function, must be sufficiently small so that the Monte Carlo statistical error is of < 1 MeV. 
A major focus of this work is on developing six CQ wave functions with small variance of 



the local energy. The exact eigenstate satisfies H^q(R.) = E ^f (R) at all R, and hence has 
zero variance. The Monte Carlo statistical errors in GFMC calculations are also determined 
by the variance of Ey(R,), thus it is necessary to obtain variational wave functions with 
small variance before proceeding to GFMC. 

The second main concern is the choice of the CQM Hamiltonian. Our H is basically 
a generalization of the Y-junction flux-tube (FT) model [OJ. The model Hamiltonian 
in Ref.|| (herein referred to as CKP) contains relativistic CQ kinetic energies, Y-junction 
FT confinement potential and one-gluon exchange Coulomb, cr^ • crj, tensor and spin-orbit 
interactions. It is not suitable to study six-quark systems for at least two reasons. 

Any model Hamiltonian used for six light quarks must be able to describe two free nucle- 
ons. In order to do that it must give E(p) = rriN+p 2 / (2m n) for the energy of a nucleon with 
momentum p <C m^. This requirement is easily acheived in non-relativistic Hamiltonians 
by choosing the mass itiq of the CQ as mjy/3. In the Hamiltonians containing relativistic 



kinetic energies it is necessary to include boost corrections |fLq , |17 | to all interactions to 
satisfy this requirement. These are absent in the CKP Hamiltonian. In order to avoid them 
we begin with the non-relativistic model. However, the CQ momenta are often greater than 
tojv/3, and it will be necessary to work with the relativistic Hamiltonian to obtain more 
reliable results. The present calculation is just the first step. 

Second, it is well known that nucleons a few fermi apart interact via the OPEP. It has 
been shown that the tensor force in the OPEP largely determines the toroidal structure of the 
deuteron [Tj| . The model Hamiltonian of CKP does not contain pion-exchange interactions. 
Within that approach, the emission and absorption of pions by nucleons may be attributed 
to breaking of flux-tubes [|19|, an effect absent in our present work. The simplest way to 
include such processes in our CQM, which we adopt, is by coupling pion fields to the CQ. 

The advantages of replacing the one-gluon exchange (OGE) interaction in CQM by OPE, 
as well as the problems associated with this have been recently discussed by Glozman and 
Riska p0| , |2~T| and by Isgur J22|. In our model we include both OGE and OPE interactions 
between the CQ. They respectively respresent the effect of virtual gluons and color-singlet 
qq pairs missing in the CQM Hamiltonian, and both should be included. A more realistic 



CQM Hamiltonian may contain additional meson-exchange terms, as discussed in Sec.[TV. 

The details of the CQM Hamiltonian used in this work are given in Sec.|I|, and the varia- 
tional wave function and VMC calculations are described in Sec.|TT|. The results, indicating 
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the limitations of the present Hamiltonian, are presented and discussed in the last section, 
Sec.[IV]. Continuing research and conclusions are given in the last SecjV]. 



II. FLUX-TUBE CONSTITUENT QUARK MODEL HAMILTONIAN 
A. Single nucleon Hamiltonian 

The Hamiltonian for the single nucleon is written as a sum of kinetic energy and two and 
three-body potential operators, 

HN = y t( m Q + Sr) + E^ + ^-Vo- (5) 

i=l ^ ®' i<j=l 

The two-body term, dij is summed over pairs (ij), and V3 is a three-body interaction. These 
potentials are due to FT confinement, and perturbative OGE and OPE between constituent 
quarks. The Vq includes all constants in the long range confining interaction and is taken 
as a parameter to fit the mass of the nucleon. In the case of the single nucleon in its 
rest frame, the constituent quark mass, rriQ, is simply a parameter which may be arbitrarily 
chosen to obtain agreement with experimental data. However, in the present non-relativistic 
QQ Hamiltonian the constituent quark mass must be one-third of the nucleon mass for the 
Hamiltonian to reduce to that of two free nucleons when two three-quark clusters are far 
apart. For this reason we use mq = in this work. 

We use the SU(3) FT model of long-range confinement The energy of the flux tubes 
connecting the three quarks is given by, 

3 

V c (r 1 ,r 2 ,T 3 ) = y/^^2r iy , (6) 

i=l 

where r^y is the length of the vector r^y = r~ ry, ry is the location of the "Y-j unction" of the 
three flux-tubes determined by the minimization of V c with respect to ry. The renormalized 
string tension, y/a, is fixed from the observed single-hadron Regge trajectories. Following 



Johnson and Thorn [^j we assume, 

E) = 2n^J , (7) 

where Ej is the mass of the state with angular momentum J. The fit obtained for the nucleon 
trajectory is shown in Fig. ([!]), The slope gives \fa = 0.88 GeV/fm; a similar analysis of A 



resonances produces essentially the same result. This value of sfo y used here is smaller than 
the 1 GeV/fm used by CKP. 

Considered as quantum strings, the flux-tubes undergo zero-point oscillations which give 
rise to a constant term, which is included in the Vo, and terms which depend inversely on 
the lengths r^y [23 1. The latter are approximately merged into the color Coulomb potential. 



We may break-up V c into two- and three-body terms, 



V c {r 1 ,T 2 ,r 3 ) = ^Y i r ij + V 3 , (8) 

i<j=l 

by defining the three-body potential: 

3 
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E^-^E^ • ( 9 ) 
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i<j=i 



This break-up is useful because the ratio of the three-body to the two-body term is < 0.154 
for all configurations, saturating when the quarks lie on the vertices of an equilateral triangle. 
The V3 term is zero when the quarks are in a line, and approaches zero when one of the 
quarks is far from the other two. 

The FT confinement model is suggested by numerical studies of lattice QCD which show 
the interaction between static quarks to be linear in the quark separation at large distances 
and independent of their spin and isospin [p5fl. The mechanism which gives rise to a linear 



confining interaction was originally studied by Wilson |26j and others using path-integral 
methods. Both lattice QCD and path-integral methods indicate that at strong-coupling or, 
equivalently at length scales of order of the characteristic length scale of QCD, believed to 



be ~ 0.5 fm |27|, the gluonic flux coalesces into tube or string-like configurations due to an 
attractive self-interaction of the gluons. 

We also make the assumption that the flux-tubes move adiabatically, the quarks remain- 
ing on their lowest-energy surface. The quarks therefore determine the locations of the 
flux-tubes and neither the flux nor the Y-junction are free dynamical variables. This seems 
to be a rather good approximation since, as shown in Ref.|28|, the lowest lying hybrid states, 



which take into account dynamics of the flux-tubes, appear to be ~ 1 GeV above the nu- 
cleon. FT topologies which have more than one Y-junction or loops of flux correspond to 



higher energy surfaces |24| and are ignored. 

As the distance between sources gets very large flux-tubes may break with concomitant 
quark- ant iquark pair creation. This process has been used to describe the two-pion decays 



of mesons || as well as decays of baryons |19| . The FT connected to a given quark % may 
break and the broken piece can reattach to the tube of a different quark j. The effects of 
such interactions between quarks % and j are to be included in the meson exchange terms in 
CQM Hamiltonian. Of these we consider only the OPE in this work. The virtual qq creation 
also renormalizes the string tension to its physical value j23 . 



The total two-quark potential Vij is given by the sum of confinement, perturbative OGE 
and OPE potentials denoted by t)?- and v^: 

^ 3 = \^r l3 +v% + v1 3 . (10) 



The OGE potential was first used in a CQM by DeRujula, Georgi, and Glashow []29j who 
showed that the splitting within single-hadron multiplets are qualitatively consistent with 
those expected from OGE interaction. The OGE potential is the QCD analogue of the 
Fermi-Breit interaction in QED. 

For point particles the OGE spin interactions are singular in nonrelativistic theory. Tak- 
ing into account the finite size of CQ, we may regulate the potentials in order to obtain 
solutions of the Schrodinger equation. We use "monopole" form factors at each vertex 

where q is the three-momentum transfer and A is the momentum space cutoff, which we 
take to be 5 fm . The color-charge density of quarks is Yukawa-like with this form factor. 
The vfj is approximated by: 

«&= E *w- T i ( 12 ) 

p=c,cr,t,l!s 

where the p = c term is the color-Coulomb interaction with 0\j = 1; p = a term denotes 
the color-magnetic contact interaction with Ofj = <Ji ■ <Tj \ t is for the tensor interaction, 
Ojj = Sij = 3crj • YijCTj ■ iij — cr i ■ <Jj\ and is is the spin-orbit contribution, Ofj = • S^. 
The operator Tj ■ Tj is the quadratic Casimir of 577(3) acting on color indices of the quarks. 
The color wave function of three quark states is antisymmetric. It is an eigenstate of Tj ■ Tj 
with eigenvalue —2/3. The specific forms of the monopole regulated potential functions: 
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are plotted in Fig. ([2]). The functions Ya and Ta are the Yukawa and tensor functions, 

— Ar 

Ya(t) = — (17) 
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(18) 



The strength of the OGE potentials is given by the perturbative strong coupling constant, 
a s . We fix the value of a s = 0.61, consistent with values in earlier works, to reproduce the 
N-A splitting. 

The OGE tensor and spin-orbit terms are not as well-established as the color- Coulomb 
and the short-ranged spin-spin terms; the latter gives a large contribution to the splitting 
of the pseudoscalar mesons, vr — p, and of the S-wave baryons, N — A, while the Coulomb 



term is seen in lattice QCD treatments of the static QQ potential |25[ . 

There is much discussion in the literature concerning the apparent lack of the spin- 
orbit term |3|, [3(1 |3l| in nucleon spectra. Many workers simply discard the longe-range, 
spin-dependent parts of OGE-the tensor and spin-orbit forces-arguing that those terms 
may be significantly modified by FT formation and are not well-supported by experimental 
data. Isgur [|I| has argued that the OGE spin-orbit term exists, but it primarily cancels the 
neglected spin-orbit interaction due to the Thomas precession in the confining FT potential. 
Oka and Yazaki |TT[ have shown that the tensor part of v\a also gives a small contribution 
to the baryon-baryon interactions. Neither the spin-orbit nor tensor parts of vfj seem to be 
important in determining the AW interaction in the state with deuteron quantum numbers. 

The OPE interaction has the form, 
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(19) 



For point particles, v™. 



SR/ 



is the 5-function part, while v^' T LR (rij) is the Yukawa function 



part of the cr^ • <TjTi ■ Tj interaction. The quark form factors modify the radial functions, 
plotted in Fig.(||), as follows: 
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<(r) = %^ 



X 



1A ' A ' (22) 



r»- ? T A (r)-~l ? -l)(Ar + l)y A (r 



Here \x = 138 MeV, the average mass of the three charge states of the pion, and and 
are as given in Eqs.( |T7| , |l8"D with A — > fi. Note that v*f R and v£' T LR have equal and opposite 
volume integrals. 

At large distances the sum of OPE interactions between the nine pairs of quarks in 

different nucleons must be equivalent to the AW OPEP tail. Thus we fix the coupling 

constant, /^qq/47t, by considering a configuration of quarks representing two well separated 

nucleons. In this limit the ratio of the sum of OPE interactions between the quarks of 

different nucleons and the v^ N contains the SU {2) spin ®SU (2) fi avor factor of 9/25. Therefore: 

3 

fnQQ = -^fwNN, (23) 

where the pion-nucleon coupling constant f^NN is known from AW scattering experiments 
to be f% NN /4:7r = 0.075|32[. Since the quoted value of the tcNN coupling constant is extrap- 
olated to pion-pole, the monopole form factor: 

A 2 - u 2 

F(q2) = WT^' (24) 

containing the pion mass \x is used to calculate the v n . We assume that the cutoff parameter 
A is the same for OGE and OPE interactions between quarks, though they could be different. 

Presumably we are justified in retaining the long-range parts, v^ LR (r) and v^ T (r), of 
OPEP. The short-range v^ T SR (r), term however, is the smeared delta-function contact po- 
tential whose validity is questionable. Within the meson exchange picture one expects an 
infinite number of qq states to contribute at short distances. The contribution of higher 
mass meson states can modify the strength of the contact term. 

The single N and A properties calculated with this Hamiltonian are listed in the rightmost 
columns of Table |I|. We note several features of these results. The kinetic energy of the 
quarks in the A-resonance is significantly smaller than that in the nucleon, thus the N — A 
mass difference cannot be calculated perturbatively as was noted by CKP. 

The v^' r SR term makes a significant contribution to the nucleon energy and the N — A 
mass difference. The nucleon is composed primarily of T = 0, S = and T = 1, S = 1 pairs, 



with nearly half the pairs in each channel. The value of <Tj • crjTi ■ Tj in these pairs are 
+9 and +1, respectively. The A however has T — 1,S — 1 pairs almost exclusively. This 
results in a strong attractive short-range pion interaction in the nucleon relative to the A 
and causes the d (u) quark in the proton (neutron) to have a smaller rms radius than the 
average. This effect may be seen in the neutron charge density, p\ , plotted along with the 
quark density in Fig.(g). 

The other short-ranged interaction is the v 9 a term in vf-. Its contributions to the energy 
of the nucleon and the iV — A mass difference are of the same order as those of vj- R . In 
fact, one can reproduce the N — A mass difference with only the v% using a larger a s and/or 
A. The fnQQ is fixed from the observed J^nn, and it is necessary to have A ~ 5 fm -1 to 
get the conventional value of a s and reproduce the A^ — A mass difference. The long-range 
spin-dependent terms, the tensor, tr, and spin-orbit, give ^ 1% contribution to the total 
energy of the nucleon. We have verified that the spectrum of P-wave nucleon and A states 
obtained with the present Hamiltonian is as good as that obtained with the semi-relativistic 
Hamiltonian excluding OPE interaction used in CKP. 

The rms quark radius for the nucleon is 0.44 fm. Taking into account the size of the CQ 
we obtain for the proton and neutron rms charge radii 0.67 fm and —0.015 fm, respectively. 
They are smaller than the observed values of 0.79 fm and —0.34 fm. A part of this difference 
could be due to the contribution of the pion cloud to the charge radius. 



B. Hamiltonian for six-quark two-nucleon states 

Possible FT configurations for six-quark states consistent with gauge invariance are shown 
in Figs.(|5]a) and (||b). The "exotic" hadron configuration of Fig.(|5|b) has been shown to lie 
~ 300 MeV above the two nucleon state and is not considered in this work. 

Fig.([|a) shows one of ten possible FT configurations with two Y-junctions, corresponding 
to one of the ten ways to divide the six quarks into two indistinguishable color-singlets listed 
in Table pTT[ Each configuration V = {1, 2, . . . , 10} specifies a color state \V) in which the 
colors of quarks ijk and Iran are coupled to separate singlets. We restrict the Hilbert space 
of the six-quark states to those of the type: 

10 

^ = J2^v\V). (25) 
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Where the ty-p are functions of the positions, spins and isospins of the six quarks. The 
ty-p\V) is antisymmetric under exchange of any pair within the singlet ijk and within Imn. 
It is also made antisymmetric under the triple exchange of ijk with Imn. The terms in the 
sum over V may all be reached from the wave function in a given partition, say V&ijl), by 
quark exchange operators, Pij shown in Table JTJ. It is useful to write Eqn.(^) explicitly 
as, 

*=(!-£ E ^pil 1 )- ( 26 ) 

\ 1=1,2,3 j=4,5,6 / 

The minus sign in the second term ensures that this wave function is completely antisym- 
metric. 

Within this space, our Hamiltonian is a generalization of that for single nucleon discussed 
above. It is given by: 

« = t ( m « - + £ («»,,+£.-„,) 

9=1 V Q / q<q'<6 

+ - (V c (Ti,Tj,r k ) + V c (r l ,r m ,r n ) 

+ V c (r t ,, tj,, r ¥ ) + V c (r v , r m /, r„0) - 2V Q . (27) 

The quarks ijk and Imn are in separate color-singlets in the ket, \^), while i'j'k' and Vm'n' 
are singlets in the bra (\&|. 

The above Hamiltonian has "P-diagonal" or color- diagonal (CD) matrix elements, 

H vv = {m v ;V\H\$ V ]V), (28) 

in which the flux tubes remain unchanged. In color-nondiagonal (CND) elements: 

H vv , = (typ; V\H\$ v r, V'), (29) 

having V ^ V, there is an exchange of the flux tubes. In QCD, when quarks in different 
nucleons are close enough they may exchange flux tubes. This effect arises due to magnetic 
terms in the lattice Hamiltonian of QCD which can alter FT paths and result in configuration 
mixing. The terms which mix configurations are proportional to the inverse of the strong 
coupling constant and, therefore, configuration mixing should be small at hadronic scales 
where the strong coupling constant is large. 

We may illustrate this point by first considering the case of infinite strong-coupling. In 
this limit, for a given configuration of quarks, the two possible FT arrangements indicated 

11 



by the sets of solid and dashed lines in Fig.(^), are orthogonal. This can be seen as follows. 
In the Hamiltonian formulation of lattice QCD the gauge sector of the theory may be written 
as a sum of two terms. An "electric" term, proportional to g 2 , the square of the strong- 
coupling constant, and a magnetic term, which goes like g~ A . The electric term simply 
counts units of flux along links of the lattice and the magnetic term may change its path 
in a gauge invariant manner. In the infinite strong-coupling limit, g — > oo, only the electric 
term remains. The eigenstates of the strong-coupling Hamiltonian are then single lines of 
flux along arbitrary lattice paths which originate on quarks and terminate on Y-j unctions 
or anti-quarks. Thus different FT configurations are expected to be orthogonal by the 



hermiticity of the strongly-coupled Hamiltonian ||34|| . Due to this orthogonality the CND 
elements of the strong-coupling Hamiltonian as well as those of the normalization will be 
zero. 

When flux sources are close the color-magnetic terms cause them to rearrange, resulting 
in FT configuration mixing. This gives rise to CND or quark exchange matrix elements, 
wherein quarks are exchanged between nucleons. A simple way to model the supression of 
mixing of significantly different FT configurations is to insert a factor into the CND matrix 
elements which falls off as the distance between the exchanged quarks increases. We have 
chosen the form: 

C(ru) = e" 7 ^ (30) 

where 7^ T is the range over which flux-tubes may be exchanged, and r« is the distance 
between the exchanged quarks i and I. More complex parameterizations for this factor, 
where ( depends on the positions of all six quarks, are discussed in Ref.[[J5]. The limit, 
7ft - * oo, corresponds to zero-range, i.e. no, quark exchange and is therefore referred to as 
the strong-coupling limit. The above ( factor is included in all CND matrix elements. We 
have considered three cases: the limit jft oo, finite jpr = 2 fm -1 , and •jft = 0. The 
jft = 2 fm -1 value supresses quark exchange for distances larger than 0.5 fm, which is of 
the order of the rms quark radius of the nucleon. 

The six-quark wave functions contain the color factors \V), however, it is simple to 
explicitly do the color algebra and suppress them in QMC calculations. In the case of three- 
quark hadrons this corresponds to replacing all the Tj ■ Tj operators by —2/3 and the unit 
operator in color space, 1q by one. In the six-quark case the color matrix elements factorize 
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from the rest: 

(V\0(R)O c \V) = ^2{V'\O c \V) f(iRe-^ |ri - r,|2 fJ„(R)0(R)f P (R), (31) 
-pi -p J 

where G(R) is a spin-isospin-spatial one or two-body operator and O c may be lc or T q -T q '. 
For V = V there is no quark exchange and r-j — = 0. When V ^ V the exchanged quarks 
i and I are uniquely determined by V and V. 

The P-diagonal color matrix elements are trivial: (V\tc\V) = 1 and T q ■ T q i = —2/3 if 
the quarks q and q' are in the same singlet, else T q ■ T q i = 0. 

The CND overlap factors 

{P\t C \P) = \Cv',v (32) 

where the phases C-p',v = ±1 are determined for the partitioning in Table |l| using the 
following rules. The overlap of \V = 1) with any other state \V' — 2, ... , 10) is given by, 

(^|l) = (l|P«|l) = +i; (33) 

independent of i = 1, 2, 3 and I = 4,5, 6. The factor | is the result for the general case of 
SU(N C ) which is N^ 1 with N c = 3. This leads to a suppression of exchange matrix elements 
in comparison to colorless objects. The remaining C-pi •p may be calculated by noting that 
a general overlap (V'\V), with V ^ 1 and V ^ 1, may be written as 

(V\V) = (l|Pi,P iTO |l), (34) 

where % and j are < 3 and I and m are > 4. In the case, i ^ j and / ^ m, 

(l|P«P im |l) = +i; (35) 



otherwise, when either i = j or I = m, 



{l\P a P jm \l) = ~ (36) 



The last case, z = j and / = m is a diagonal element. 
The CND (P' ^ P) factors are: 

,(2) 



cZ,p>,p = (^IVW (37) 

= (i'/jf; l'm'n'\T q , ■ J q \ijk; Iran) (38) 
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We note that pl p = Cff q . pp , = Cf^.p. p. Using cyclicity: i'j'k' = j'k'i' = k'i'j' and 
Vm'n' = m'n'V = n'l'mf we can always bring the above matrix element into the form: 

Cf q . pip = (lj"k"; im"n"\T q > ■ T q \ijk; Imn) (39) 

There are four possibilities in Eq. (|39| ) corresponding to j" = j or k and m" = m or n which 
can differ from 

{Ijk; imn\T q i ■ T q \ijk; Imn) (40) 

only by a phase given by: 

{lj"k"; im"n"\T q i ■ T q \ijk; Imn) = (2(5j»j — 1) (2<5 m » im — 1) (Ijk; imn\T q r ■ T q \ijk; Imn). (41) 

This allows us to work with just the color factor in Eq. (flCf) . There are four types of quark 
pairs q',q: (1) For the exchanged pair q',q = i,l (V'\T q i ■ T q \V) = |. (2) For unexchanged 
pairs in the same singlets, qq' — jk and mn (V\T q > ■ T q \V) — — |. (3) For unexchanged 
pairs in different singlets, qq' = jm,jn, km and kn (V'\T q i ■ T q \V) = |. (4) The remaining 
eight pairs are between an exchanged quark q = i or I and one of the unexchanged quarks 
q' = j, k, m or n. For these (V\T q > ■ T q \V) = — §. 

III. VARIATIONAL WAVE FUNCTION 

We use the variational method of Ref.[|36], applied there to light nuclei, to the problem 
of six interacting quarks. The idea behind this method is that interaction terms in the 
Hamiltonian will induce correlations in the wave function which, in general, depend on the 
quantum numbers of the interacting quarks. A good approximation to the ground state 
eigenfunction may be obtained by applying two and three-body correlation operators with 
the same operator structure as the interaction terms in the Hamiltonian to an uncorrelated 
wave function. 

The six-quark wave function corresponding to two uncorrelated nucleons is just a product 
of two single-nucleon three-quark states. We assume that the full correlated wave function 
is obtained by operating on it with correlation operators for pairs of quarks in different 
nucleons, as well as for the nucleon pair. 
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A. Single nucleon wave function 



The uncorrelated nucleon wave function, § N (m,T,ms) is the product of symmetric spin- 
isospin state with T — S — |, T z — uit, S z = m s , and antisymmetric color-singlet wave 
functions. The latter is eliminated from QMC calculations as described in part (B) of the 
last section, and the former is, for a spin-up proton, for example: 

$iV (+7T> = ~7= ( 2 I M J [u J [dl)-\u J [uid J [)-\uiu J [d J [) + perms.) . (42) 



2' 2) 3V2 

It has no dependence on the quark positions. The variational nucleon wave function used 
in this work is, 

^ N (m T ,m s ) = F/ 23 ^n^) ® N (m T ,m s ). (43) 

Here, F$ is a three-body correlation function and are pair correlation operators. The 
superscript / distinguishes correlations between quarks internal to the nucleon. When we 
consider six-quark wave function, it will have another correlation operator, F?, which acts 
on quarks in different nucleons. The symmetrized product is required since F^ do not 
commute. The ^f N is translationally invariant. 

In CKP the CQ pair interaction contains only the confining term and the OGE interac- 
tion. There the form of F/, is taken to be, 

F\r P) = I 1 + Mrij)<Ti ■ ctjUMj), (44) 

where f c denote spatial correlations and u a (rij) is the spin-spin correlation function induced 
by the <Tj • crj term in the OGEP. The correlations induced by the tensor and spin-orbit 
parts of the OGE interaction are small, and neglected in CKP. 

In the present work, the CQ interaction also includes the OPEP. We therefore consider 
the most general static spin-isospin correlation operator: 

6 



p=2 



fcira), (45) 



where the sum runs over the operator designations: r, a, err, t, and tr which correspond to 
the operators, numbered p = 2 — 6, 



Oij 2 6 = Tj • Tj,<Ti ■ <7 j,<7{ ■ (TjTi ■ Tj, Sij, SijT j • Tj . (46) 
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The Ojj = 1 is denoted by symbol c. The r term is included for completeness, even though 
it does not appear in the OGE and OPE interactions. We neglect the spin-orbit correlations 
since the spin-orbit interactions play a small role in the present problem. 

The three-body correlation function takes into account the small, spin-independent three- 
body potential of Eq.(^). We take the functional form of the correlation suggested by first 
order perturbation theory, 

^23 = V3(ri,r 2 ,r 3 ) , (47) 



where f3% is a positive variational parameter. In fact, we find that the value used in Ref.p5 
A, = 0.025 x 10~ 3 MeV" 1 is sufficient. 

The pair correlation functions, f* (r#) and uUrij) = fp{rij) / f* (r#) are varied to minimize 
the energy of the single nucleon states. Since the pair interaction (Eq. [10]) depends upon the 
total isospin and spin of the interacting quarks it is convenient to project the correlations 
into the four possible T, S channels, 

4 

AsMlT, S) = J2 fti^O^T, S) , (48) 
P =i 

flAr^T, S)=J2 fp( r v)°ii\ T > S ) ■ ( 49 ) 

p=5,6 

In a single nucleon, the antisymmetry is ensured by the color-singlet factor, implying that 
the spin-isospin spatial part of the wave function is symmetric. Thus the two-quark states 
have T,S — (0, 0) or (1, 1), for the smallest L = 0, while L = 1 for T, S = (1, 0) or (0, 1). 

Let vip $ an d v\ T be the projections of the pair interaction (Eq. [10]). In absence of the 
third quark the pair correlation functions would obey the following two-body Schrodinger 
equations with a constant A^ 5 , representing the eigenvalue, and \ tT = 0. The equations in 
the S = channel are 

- - ^P^ir^fT,*)) + [4,o ~ A^ ](r L+1 /^o) = 0, (50) 

with L — 0, 1 for T — 0, 1 respectively. The coupled equations for S = 1 are, 

+v / 8K 1 -A( 1 ](r^r) = (51) 

-^[(r L+1 fl T )" - *- ± ^ L {r L+1 fl T )] + - 4,i - 2«t - Kt) ~ Hs,T)(r L+1 fl T ) 

+V8[vl T -\U(r L +ifi)=0 (52) 
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with L = 0, 1 for T = 1,0. The L,S = 0,0; 1,0 and 0,1 equations are for 1 5'o, l P\ and 
3 Si — 3 D\ waves, while that for the L, S — 1, 1 is an average for the 3 Pj states as discussed 



by Wiringa [5c 



The effect of a third quark on the interacting pair is taken into account by taking the X 1 
as functions parameterized as: 



= + (l - e-<*,>°) (53) 

a j,o 

\! (r) = , . (54) 

The A 7,0 's are adjusted to match the boundary conditions discussed below, and the second 
term in Eq.(|53[) is required to obtain the boundary conditions. The parameters, Rt,s, Rt,Ti 
aT,s, <H,Ti an d £,t,s are varied to minimize the single-nucleon energy. We reduced this set 
of 16 variational parameters to six, assuming Rr,s = Rt,T = Rs, a r,s = a>t,T = «5 and 
£t,s = £s- The small statistical variation of the total energy of the single-nucleon state 
suggests that this set of six parameters is sufficient. 

The equations are solved with boundary conditions appropriate to the three-body system. 
When one quark is far from the other two, the linear confining potential dominates in the 
three-body Schrodinger equation which reduces to, 

lim i^V 2 (/f(r)) 2 - vM/f(r)) 2 ) = 0. (55) 



The solution is the asymptotic form of the Airy function, r (//(r)) 2 — > e Kri/2 with k = 
\6my/o : /27h 2 . The full set of boundary conditions is then, 

/t s( r 0) = constant (56) 

#, s (r-oo) = /^W— — (57) 
flAr - 0) = (58) 
fl T (r - oo) = i*r/£i(r) (59) 

where the central normalizations, /?-t,5, and D/S ratios, T]t, are further variational param- 
eters. Only three normalization factors, hT,s, are independent since one sets the overall 
normalization of the wave function. 

The optimum internal correlations are plotted in Fig.(|7]). At small distances, < 0.5 fm the 
spin-isospin dependence is seen in the relative sizes of the s . We notice that T, S = 0, 
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pairs are favored at small distance due to the attractive short-range terms in the OGE 
and OPE interactions. This results in the w-quark (d-quark) distributions in the neutron 
(proton) having smaller radii than their isospin partner and gives a negative rms charge 
radius to the neutron, as mentioned earlier. 



B. Two-Nucleon variational wave function 



The six-quark wave function representing two uncorrelated nucleons in T, S = 0, 1 state 
with spin projection Ms is given by: 

10 

$ NN (M S ) = ^(-l) (1+a ''- l) $£ JV (M s )|7>) (60) 
v=i 

$% N {M S ) = -^[^(p,m 1 ;ijk)^{n,m 2 ;lmn) 
v2 

- ^p(n,m 1 ;ijk)^p(p,m 2 ;lmn)] S M '^ (61) 

Here ^f^(p(n),mi;ijk) denotes the wave function of quarks ijk in the proton (neutron) 
with spin projection mi, as in Eq . (|43|) . The quark numbers ijk and Imn depend upon the 
partition V, and the full wave function Q NN (Ms) is antisymmetric under the exchange of 
any quark pair. The spins inside square brackets are coupled to values indicated by the sub- 
and superscripts. 

The variational six-quark wave function, \lV for the interacting two-nucleon state is 
obtained by inserting correlations in the above wave function. We consider long range 
correlations between the centers-of-mass of quarks ijk and Imn, as well as correlations 
between the nine external quark pairs, like il, with one quark from each singlet. The 
correlations between internal quark pairs, such as ij in the same singlet, are included in the 
nucleon wave functions ^p(p(n),mi;ijk). The is chosen as: 

/ \ 



m v (M) = ^(-i)(!+^.i) 



V=l 



s n *5 



q = i,j,k 
q 1 =1 ,m,n 



*$ N (M)\V) , (62) 



9™(M) = ^Y ^ N (Ms = M) + _i_L-Z. [Y m, (? V )*»»(M 8 )] "~_ M (63) 

T-p T-p J 



where 

i> = o( r « + r i + r k - ri - r m - r n ). (64) 
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The long range correlations, u(r-p) and w(r-p), are calculated from the two-nucleon 
Schrodinger equation, 

■— V 2 + v NN (v) } * NN (v) = E NN * NN (r). (65) 
m N J 

Here V NN is the wave function of Eq.fl6"3|) taken in a given partition, say the first. The 
substitution of ty NN into Eq.(^) results in the set of coupled equations for u and w shown 
in Eqs.([I],|2|) with u — > it and w — > w. These are solved subject to the boundary conditions: 

u(R c ) = 0. 

w(R c ) = 0. (66) 

for the ground state in a cavity of diameter R c . The two-nucleon effective potential v^n, 
used to minimize the energy, is discussed in the next section. 

The external correlation operators, F^,, contain six terms associated with the operators 
O^ 1 ' 6 , as do the the internal correlations. By projecting the F^ q , into two-quark T, S 
channels we obtain the six functions f^ s and ff T . These functions are assumed to obey 
Eqs. ([50H52|) with the internal v 1 and X 1 replaced by the external v E and X E , however the 
boundary conditions of the f E are different from those of f 1 as discussed below. 

Terms without quark exchange between the nucleons, i.e. those having V = V seem to 
dominate the expectation values of 1 and H. In these terms the confinement as well as the 
OGE interactions do not contribute to external pairs. Therefore we take: 

The boundary conditions for the external correlations are those appropriate for objects 
which are free at large separation. They are: 

f% 8 (r -> 0) = constant (68) 

fis(r > 4) = 1 (69) 

f t E T (r - 0) = (70) 

f t E T (r > dt) = (71) 

where d c and dt are the central and tensor "healing" lengths which are varied to minimize 
the variational energy. 
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As in the theory of nuclear matter |37[ the X E are taken as constants for distances less 
than the healing lengths while the boundary conditions fix the form for radii greater than 
the healing lengths: 

A|, (r) = Xrp(dc -r) + v Tfi 9(r - d c ) (72) 
Af,i(r) = 4$6{d e - r) + [v$ tl + V8(v E T - X E A )f t , T ]9(r - d c ) (73) 
>% T (r) = X?ft(d t -r) + v t>T e(r - d t ), (74) 

with 6(r) the Heaviside step function. The X E, ° are varied to match the boundary conditions, 
so there are no free parameters besides d c and d t . The external correlations are shown in 
Fig.(|) for the optimum values of d c , d t at R c = 2 fm. The central correlation, f E {r) is seen 
to be ~ 1 over the entire range, with a slight 'dip' at short distances which arises primarily 
due to the short-ranged err term in OPEP. Even though the f E are small, they have an 
effect since the total wave function has a symmeterized product of nine of them. 



C. Monte Carlo evaluation 



The present VMC calculation of six-quark states is formally similar to that of 6 Li [13 
except for the color factors. The wave function Eq . (|62"D is expressed as a vector function of 
the configuration vector R = {r 1; ...r 6 }: 
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10 320 

= EE^v*( R )i a >i p >' ( 75 ) 

V=l a=l 

where \a) are six particle spin-isospin states. The tensor correlations populate all the 2 6 
spin states, and, since isospin is conserved, we consider only the five possible T = six- 
quark states. The set of the 320 states |a) used here is the same as that in 6 Li. In this 
representation the correlation operators F qq i are sparse, 320 x 320 matrix functions of R. 
The wave function of Eq.(|62]) contains symmetrized products of correlation operators F^ q , 
and F q E q/ . Each * P (R) is a sum of A^ r ! terms, where the number of pairs, N pr = 15. The 
^/•p jr denotes the term in which the N pr pair correlation operators act in an order denoted 
by r, and 

N pr \ 

vMR) = J^, r (R). (76) 



r=l 
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This sum is evaluated stochastically in Monte Carlo. 
The expectation value of an operator: 





0\V V ) 




$v) 



is evaluated using N c samples of R/, and left, right operator orders lj and rj obtained by 
sampling the weight: 

Wi, r (R) = | ^(-l)^- 1+ ^M^(R)^ ir (R)(P|P / )e~^|r l -r J f| ) ( 78 ) 

v,v 

The average value and the standard deviation are calculated from ~ 10 5 samples using 
block averaging. The CD and CND contributions are evaluated separately; for example 
N = Mcd + Ncnd, where Men is the average value of: 

^^(R^O^R,) , (79) 
v 

while Ncnd contains all the other terms. 

The Fortran90 computer code for this calculation was run on 6 nodes of an IBM-SP . The 
amount of computer time required to generate one statistically independent configuration 
{R, l,r} and evaluate its energy is 1.35 node seconds. For runs of 10 5 configurations, we 
obtain a statistical variance for the total energy of the system at the ~ 0.1% level. A run 
of this length takes ~ 40 node hours. 



IV. SIX-BODY RESULTS 

The parameters for the single particle Hamiltonian, Eq.([|), as previously described, are 
listed in Table |V[ As mentioned earlier the a s is determined from the N — A mass difference 
and Vq is chosen to reproduce the nucleon mass, m N = 938.9 MeV, taken as the average of 
the proton and neutron masses. We choose two values for the flux-tube overlap parameter, 
7ft, to cover the physically interesting range, and minimize the energy of \lV with respect 
to variational parameters for each of them. The energy is also calculated for the R c = 2 fm 
case, with 'Jft = 0, perturbatively, as discussed below. 

The six-body wave function is minimized with respect to the external pair correlation 
operators, F^, and the iViV long range correlation functions, u and w. The Fy are varied 
by adjusting the central and tensor healing lengths, d c and d t . Internal pair correlation 
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operators, Ft, as well as the internal triplet correlation function, F( 23 , are held at the 
equilibrium values for the single nucleon. The u and w in Eq . (163|) , are obtained from Eq . (|65|) 
with an effective NN interaction, v^n and subject to the boundary conditions, Eqs.(^). 
We obtain an approximate effective potential, v con by the "convolution" method: 

v con {R; M)= [dK S^U - |z) <^(r 2 + |z) $r f (M) ( £ *«) 

^ ' ^ ' V 9=1,2,3 g' =4,5,6 / 

(80) 

with <iR = d 3 rid 3 r2 ■ ■ ■ d 3 r$, Ri = |(ri + r 2 + r 3 ), R 2 = |(r 4 + r 5 + r 6 ), and <&i N given by 
Eq.fl6T|) with V = 1. The delta-functions fix the centers-of-mass of the nucleons at separation 
R along the z-axis. Note that <$>^ N depends only on the internal quark pairs in the first 
partition V = 1 and is independent of the nucleon centers-of-mass, Ri and R 2 . Only the 
OPEP between quarks in different nucleons contributes to this convolution potential; it 
therefore has no spin-orbit term. We evaluate v con (R; M) via Monte Carlo integration for 
two values of M = 1, and project into central and tensor potentials, v c con and v \ on as: 

< on (R) = l l2v con (R; i) + v con (R; o)] (81) 
vlniR) = I Kon(R; i) - v con (R; o)] . (82) 

These are plotted in Fig.(^) along with the corresponding terms from the Argonne fig 
potential which appear in Eqs. (|l],0) . 

The effective potential appearing in Eq. (|65|) is taken as: 

v NN = a c v c con + a t v t con S 12 - (83) 

where a c and a t are determined variationally for each R c and ^ft- The resulting u and w 
for '-/ft — 2 fm" 1 are plotted in Fig. (|TDD . They should not be interpreted as the iViV two 
nucleon relative wave functions since the six-quark \Ev contains additional f E correlations. 
The iViV radial wave functions may be obtained from tyy and ty NN , as discussed in Sec.0. 
Table ^ gives the optimum values of d c , d t , a c , and a t . As discussed later, our \Ev is less 
accurate for R c = 6 fm, jpj< — > oo than for all other cases. The a c and a t for this case may 
not be the optimum values. 

Smaller values of a c are favored presumably because the external pair correlation function, 
fc{ r ij)-> l ess than one for small distances, supresses configurations where the quarks in 
different nucleons approach each other. When the nucleons overlap the small difference in 
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fci r ij) from unity has a significant effect on the wave function since it contributes via all 
the nine pairs of quarks belonging to different nucleons. The facts, that at ~ 1 in all states 
and d t is small, suggest that most of the tensor two-nucleon correlation is induced by the 
effective vnn- 



top three rows of this table give the total, kinetic, and potential energies for the optimum 
\EV per nucleon. The statistical variance relative to the total energy is: 



where £y(Rj) is the local energy, Eq.(||), Ey is the average of the local energy, and N c is 
the number of configurations, R/. The relative fluctuation, e(R/), in Ey(Ri) is defined as 
Ey(Ri) = Ey{l + e(R/)), and a re i gives its rms value. If \lV is the true ground state wave 
function then EV(Rj) = Ey, for all Rj and e = 0. The rms value of e is 0.13 in all states 
in Table [il] except 'Jft - > oo, R c = 6 fm where a re \ = 0.17. The perturbation calculation, 
R c = 2 fm, jft = has o re \ = 0.19. We can probably reduce these e values by reoptimizing 
a c , act, d c , and d t , but the present results serve to indicate the role of the external f E and u 
and w in improving \&y. The fluctuations in the kinetic and total potential energies, T and 
V, relative to the total energy, are much larger, of ~1. Although T and V fluctuate wildly 
as R/ varies, their sum, Ey(Hi), changes little, evidence that \Ev is a good representation 
of the true ground state, ^o- 

The standard deviation in the expectation values obtained from 10 5 configurations is 
shown in parentheses for all observables in Table [Tl], and the deviation in Ey is only 0.4 
MeV per nucleon. This table also lists the contributions of the various potential terms in 
the Hamiltonian, and the energy of \Ev relative to two free nucleons, per nucleon, compared 
to the empirical value calculated from Argonne t>ig. It shows that all of the terms in the 
Hamiltonian differ by only a few percent from their single nucleon values, except v^ T , which 
is greatly enhanced. Comparison of the calculated 5E(R C ) with 5E emp (R c ) shows that the 
present six-quark model does not give a sufficiently attractive strong interaction. The quark 
kinetic energy gives a repulsive contribution to E(R C ), while the interaction terms give a 
much smaller attraction. For example, at R c = 2 fm the iViV values for the total kinetic 
and potential energies per nucleon in Table | are 80 and —43 MeV, respectively, while the 
six-quark values are 62 and —3, relative to two free nucleons. The kinetic energy is of the 



The results of the full VMC calculation are given in Table [TT] for each R c and 7ft- The 




(84) 



23 



expected size but the potential is not attractive enough. 

Comparison of jft = 2 fm -1 and oo results indicates that the total effect of CND contri- 
butions, via flux-tube or quark exchange, forbidden for jpr — > 00, is small and repulsive. A 
perturbative calculation with the jft — 2 fm -1 wave function for jft = 0, i.e. no exchange 
suppression due to flux-tube overlap, gives an even higher energy of 1015.0(6) MeV/nucleon 
at R c = 2 fm, corresponding to E(R C ) = 76.1(7) MeV/nucleon and 5E(R C ) = 24.9(7) 
MeV/nucleon. We breakdown the expectation values of operators into CD and CND terms, 
listed in Tables |V H[VI]| and [VI I H , according to: 

(O) = ^"If"" (85) 

L^CD + UCND 

where the terms in the numerator and denominator are evaluated as in Sec. [III C[ The CND 
term of the unit operator, Vcnd, is small so we expand the denominator of Eq.(|85|) and 
obtain the CD contribution as Mcd/T^cd-, while the CND part of the expectation value is: 



NcND T^CND A/"cd 



(86) 



T^CD T^cd T^cd 

As can be seen from Tables |VT| to |V111| , the variance of the CD and CND contributions to 
the energy is larger than that of the total. 

For short range quark exchange, ^ft — 2 fm -1 , the CND contribution of the unit operator, 
T-^cnD) is —0.0142, —0.0017, and —0.0003 with the wave function normalized to one, for 
R c = 2, 4, and 6 fm, respectively. However, the sign of the T>cnd changes for long range 
quark exchange, jft — 0, R c = 2 fm, having the value +0.0058. It is simple to understand 
why the exchange contribution is expected to be small and positive for long range quark 
exchange, so we consider this case first. If we ignore all spatial correlation of quarks in 
the wave function Q NN (Ms) of Eq. (|6fj|) , we obtain a CND contribution of +1/9 to the 
normalization. Localization of the nucleons reduces this value by suppression of exchanges 
for separations of the exchanged pair larger than the single nucleon rms quark diameter. 

In the case of jft — 2 fm -1 the sign of the CND contribution to the normalization is a 
consequence of the spin-isospin-color dependence of the wave function. Overall antisymmetry 
of the wave function implies that the sign of the contribution for exchanged quark pairs in 
the T,S = (0,0) and (1,1) states is negative while those in (1,0) and (0,1) are positive. 
External quark pairs interact via OPEP, cf. Eq. (|H7|) , resulting in the strong spin-isospin 
dependence of the external pair correlation functions, f^ s , shown in Fig. (p+f ). The is 
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attractive at short range in T, S = (0, 0) and (1, 1) states resulting in the enhancement of 
exchanges which make negative contributions to the CND normalization. Similar reasoning 
shows that the exchange contribution due to T, S — (1, 0) and (0, 1) pairs is suppressed. 

Though the quark exchange contributions are < 1% of Ey, they are ~ 10% of the 
E(R C ), which is related to the two-nucleon interaction, and add to its repulsive core. They 
significantly effect the total effect of strong interactions, being the same order as 5E(R C ). 



V. CONCLUSIONS 



We have shown that it is possible to obtain a variational wave function for the six- 
quark system which gives a small standard deviation in the energy for relatively modest 
amounts of computational time. The present study also shows that exact calculations with 
small statistical errors are possible with this Hamiltonian. A GFMC code is currently 
being developed to project out the true ground state component of using Euclidean 
propagation, ty = lim r ^ 00 e~ ( - H ~ Eo ^ T ^fy. 

This study also highlights the need for an improved CQM Hamiltonian. A primary issue 
is the inclusion of relativistic effects. The relativistic kinetic energy, T = £V + m ? is 



easy to include. Boosting the two-body potentials has been studied by Carlson et. al.fl6 



and Forest et. al.||17| in the context of light nuclei, and by Isgur |38| in a CQM. The rel- 
ativization of the current model may have a significant effect on the six-quark NN state. 
One consequence of including relativistic kinetic energy, as shown in CKP, is that with it 
nucleons have a smaller rms quark radius. Using the parametrization for ty N from CKP, we 
obtain an rms radius of 0.33 fm, against the 0.44 fm of the present *f? N . This has a significant 
effect on the convolution potential, Eq.(^), as seen in Fig.(^). As the rms quark radius of 
the nucleon goes to smaller values the central term of the convolution potential gets more 
repulsive near r ~ and, more importantly, the tensor gets stronger. Fig.(^|) shows the 
convolution potential tending toward the Argonne t>ig curves as the nucleon quark radius 
decreases. 

The meson exchange picture, correct when nucleons are at large separation, is question- 
able when the nucleons overlap considerably. While the v^ T LR and v^ T are fixed by the data, 
the validity of v*f R is questionable. As seen from the tables, this term gives a large repulsive 
contribution to 5E(R C ) and one of the ways to get better agreement with 5E emp (R c ) may 
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be to reduce its strength. 

Accepting the limitations of the present model and calculations we can add terms to the 
six-quark Hamiltonian to reproduce the 5E emp (R c ). For example, if we assume quark pairs 
interact via a medium range attractive scalar potential v s of the two-pion exchange range: 



v s = c s 



q<q'<6 

^-Mr As -^(M-i)(A sr + 1) n s , 

its strength c s and cutoff As can be adjusted to fit 5E emp (R c ). As before /i is the pion mass. 
The two-pion exchange interaction between quarks has recently been discussed by Riska 
and Brown Jj9] . The effect of this potential is calculated perturbatively from the two-quark 



distribution functions, 

P2\f) — i 88 

which give the probability to find two quarks a distance r apart. They are plotted in Fig. (|T2"D 
for R c = 2, 4, and 6 fm. The expectation values of v are given by: 



(v b ) = J d i rp 2 (r)v b (r). (89) 
The results of this exploratory calculation are shown in Table for As = 5 fm" 1 and c s 



chosen to reproduce SE emp (R c = 2 fm). The (v ) in the single nucleon state is —51.4 MeV 
and —21.0 MeV, for A = 5 and 2 fm -1 cases, respectively, which is small enough to be 
treated perturbatively, and it merely redefines the constant Vq. Though the effect of v s is 
minimal on the single hadron spectrum, it is larger in the six-quark case since there are 
nine pairs of quarks in different nucleons. This illustrates how we could refine the CQM 
Hamiltonian using NN scattering data. 

The quark pair distribution function in a single nucleon, (r), is also plotted in Fig.(|T2"D, 
scaled by a factor of 2 for easier comparison with the p 2 for six quarks. When R c = 4 and 
6 fm the first peak in p 2 at r ~ 0.6 fm corresponds to the distribution of pairs in a single 
nucleon while the second peak at larger radii corresponds to the distribution of pairs of 
quarks in different nucleons. In these cases the first peak approximately equals 2p 2 (r) as 
expected. For R c = 2 fm these peaks merge. 
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Further insight into the quark substructure of the two-nucleon states is gained by calcu- 
lating the quark-pair distribution function for the T q ■ T q > operator in OGE interactions. It 
is denoted by, pf, and defined as: 



l>2 '' /dR*{,(R)* v (R) 



and plotted in Fig.(|T^). The figure has four overlapping curves showing f° r Rc = 2, 4, 6 
fm and 2p^ N ^ in the single nucleon. The T q ■ T g i has the value —2/3 for quark pairs single 
nucleons and is zero for CD terms when the quarks are in different nucleons. Considering 
only the CD terms we expect the p\ to lie on top of each other for all R c . CND contributions 
can cause small (~ 1%) differences in p\ dependent on R c . However the accuracy of the 
present calculation of p% if) is also ~ 1%. Note the absence of a second peak for r > 2 fm 
indicating the lack of long-ranged color-dependent correlations, as expected. 

Although we will not include the results here, the effective NN wave functions may be 
calculated from the tyy using the methods of Schiavilla et. al.[40[ and Forest et. al.[|18[], and 
the knowledge of E(R C ) allows one to determine NN scattering phase shifts ||41|| . The aim 
of the present work was to explore if QMC calculations can be used to study the six-quark 
system with sufficient accuracy to discern nuclear effects. The present results are obviously 
encouraging. 
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R c 


2 f m 


4 fm 


6 fm 


oo 


E(R c )/2 


37.01 


2.32 


-0.35 


-1.12 


(T N (R c ))/2 


79.61 


23.96 


15.16 


9.94 


(v NN {R c )}/2 


-42.62 


-21.64 


-15.51 


-11.06 


Eni{R c )/2 


51.16 


12.79 


5.69 





5E emp (R c )/2 


-14.16 


-10.48 


-6.04 


-1.12 



TABLE I: Energies per nucleoli calculated from the isoscalar part of Argonne vis potential (in 
MeV) vs. the diameter of the confining cavity. R c — > oo corresponds to the free deuteron. 




1/2 3/2 5/2 7/2 9/2 11/2 

Spin 



FIG. 1: Nucleon mass squared versus J in GeV 2 . The slope is reproduced with string tension 
= 0.88 GeV. 
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R c 


^FT = 2 fm 
2 f m 4 fm 


-l 

6 fm 


2 fm 


iFT ~ * 00 
4 fm 


6 fm 


N 


A 


Ey 


qqs q(a\ 


q^o o(a) 


943.3(4) 


990.7(4) 


948.1(4) 


942.6(5) 


938.9(4) 


1233.2(4) 


T 


1 1 1 r i(A\ 


1U I Ul 1 


1062(4) 


1118(4) 


1065(4) 


1057(4) 


1053(4) 


755(3) 


V 


ZO 1 (4J 


ZDo(4 J 


255(4) 


247(4) 


257(4) 


259(4) 


260(4) 


852(2) 


v c 






995(2) 


997(2) 


997(2) 


996(2) 


996(2) 


1157(2) 


g 

Vc 






-340.8(4) 


-340.6(4) 


-340.3(4) 


-340.6(4) 


-340.8(4) 


-308.4(4) 


q 
Va 


— 101.8(3) 


— 106.5(4) 


-107.5(4) 


-107.0(4) 


-106.7(4) 


-106.9(4) 


-107.4(4) 


56.6(2) 


g 




—A (~\ 


-4.23(1) 


-4.32(1) 


-4.25(1) 


-4.23(1) 


-4.24(1) 


-7.75(7) 


g 
is 


— 9 90f 1 


— 9 0Ql4 

i.UiJI 1 / 


-2.06(1) 


-2.24(1) 


-2.09(1) 


-2.07(1) 


-2.08(1) 


-2.28(5) 


V *,SR 


— ZOO[L ) 


one/ 1 "\ 
— Zyo(l ) 


-298(1) 


-292(1) 


-295(1) 


-297(1) 


-299(1) 


-37.2(2) 


n,LR 
Var 


20.91(3) 


22.52(4) 


23.13(4) 


21.51(3) 


22.54(4) 


23.10(5) 


23.69(4) 


3.96(1) 




-27.7(1) 


-14.8(2) 


-9.7(1) 


-24.6(1) 


-14.4(2) 


-9.49(2) 


-6.48(1) 


-9.86(9) 


E(R C ) 


60.0(6) 


11.1(6) 


4.4(6) 


51.8(6) 


9.2(6) 


3.7(7) 






5E(R C ) 


8.8(6) 


-1.7(6) 


-1.3(6) 


0.7(6) 


-3.6(6) 


-2.0(7) 






$E em p(R c } 


-14.21 


-10.5 


-6.0 


-14.21 


-10.5 


-6.0 







TABLE II: Energies per nucleon (in MeV) of six-quark states compared with the single hadron 
values for single N and A. Statistical errors are shown in parenthesis. 
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V 


(ijk) (Imn] 


Pin 


1 


123 


456 




2 


423 


156 


Pu 


3 


523 


416 


Pl5 


4 


623 


451 


Pl6 


5 


143 


256 


P2A 


6 


153 


426 


P25 


7 


163 


452 


i"26 


8 


124 


356 


^34 


9 


125 


436 


^35 


10 


126 


453 


-P36 



TABLE III: Six-body wave function partitions and quark exchange operators. The first column 
is the partition number, the second and third lists the quarks joined by flux tubes in singlets, and 
the last column is the quark exchange operator which gives that partition when applied to \V = 1). 



^ 880 MeV 
a s .61 
fnNN/^lT .075 

m Q 313 MeV 

A 5 fm _1 

~/FT 2,oo fm _1 

V 373.7 MeV 

TABLE IV: CQM Hamiltonian parameters 
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R c 


7ft (fm l ) 




d c (fm) 


d t (fm) 


2 


2 


.10 .90 


3.55 


1.25 




oo 


.17 .90 


3.35 


1.10 


4 


2 


.30 1.0 


10.0 


1.20 




oo 


.50 1.0 


10.0 


1.10 


6 


2 


.90 .93 


10.0 


1.20 




oo 


.34 .72 


10.0 


1.20 



TABLE V: Optimum values of the variational parameters. 





1FT 


= 2 far 


-i 


7ft ~~ ¥ °o 




7ft = 






CD 


CND 


Total 


Total 


CD 


CND 


Total 


Ey 


990.8(8) 


8.1(3) 


998.9(4) 


990.7(4) 


990.3(8) 


25.1(3) 


1015.0(6) 


T 


1117(4) 


-2(1) 


1115(4) 


1117(4) 


1113(4) 


7(1) 


1119(4) 


V 


247(4) 


10.4(2) 


257(4) 


247(4) 


252(3) 


18.1(3) 


270(3) 


v c 


997(2) - 


0.44(3) 


996(2) 


997(2) 


1000(2) 


2.1(4) 


1002(2) 


v% 


-340.5(4) 


-.2(2) 


-340.8(4) 


-340.6(4) 


-340.0(4) 


-.5(3) 


-340.4(4) 


q 
V a 


-106.2(3) 


4.4(2) 


-101.8(3) 


-107.0(4) 


-106.0(3) 


8.0(3) 


-98.0(3) 


q 
V% 


-4.31(1) 


-.43(1) 


-4.74(1) 


-4.33(1) 


-4.30(1) 


-1.40(3) 


-5.69(2) 


q 

V ls 


-2.26(1) 


.06(1) 


-2.20(1) 


-2.24(1) 


-2.25(1) 


.024(6) 


-2.22(1) 


TT,SR 

Var 


-290(1) 


8.1(5) 


-283(1) 


-292(1) 


-290(1) 


12.9(4) 


-277(1) 


TT,LR 

Var 


21.47(3) 


-.56(3) 


20.91(3) 


21.52(3) 


21.44(4) 


-1.71(6) 


19.71(4) 


Vlr 


-27.35(1) 


-.39(1) 


-27.7(1) 


-24.6(1) 


-27.3(1) 


-1.29(3) 


-28.6(1) 



TABLE VI: Variationally determined six-quark state energies for cavity diameter R c = 2 fm for 
finite (7ft = 2 fm -1 ) and strong (7ft — ► oo) coupling. The columns give color-diagonal, CD, 
color-non diagonal, CND and total contributions. The 7ft = state is evaluated perturbatively. 
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^FT = 2 fm 


-i 


7ft ~~ * °° 




CD 


CND 


Total 


Total 


Ey 


948.8(5) 


1.2(1) 


950.0(4) 


948.1(4) 


T 


1070(4) 


-.3(4) 


1070(4) 


1065(4) 


V 


252(4) 


1.5(1) 


253(4) 


257(4) 


v c 


995(2) 


0(1) 


995(2) 


997(2) 


q 
Vc 


-340.8(4) 


-.04(9) - 


-340.9(4) 


-340.3(4) 


q 


-107.1(4) 


•6(1) - 


-106.5(4) 


-106.7(4) 


g 
v t 


-4.25(1) 


-.06(1) 


-4.31(1) 


-4.25(1) 


q 

zs 


-2.10(1) 


.002(1) 


-2.09(1) 


-2.09(1) 


7T,SR 


-296(1) 


1.1(2) 


-296(1) 


-295(1) 


w,LR 
Var 


22.59(4) 


-.07(1) 


22.52(3) 


22.54(3) 


Vtr 


-14.8(2) 


-.04(1) 


-14.8(1) 


-14.4(2) 



TABLE VII: Six-quark state energies for cavity diameter R c = 4 fm. Notation identical to Table 



400 




r (fm) 



FIG. 2: Regulated OGE potential in MeV. 
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1ft 


= 2 far 


-l 


IFT 0° 




CD 


CND 


Total 


Total 


Ey 


943.1(5) 


.27(3) 


943.3(2) 


942.6(5) 


T 


1062(4) 


.0(1) 


1062(4) 


1057(4) 


V 


255(4) 


.23(3) 


255(4) 


259(4) 


V c 


995(2) 


-.01(5) 


995(2) 


996(2) 


q 

Vc 


-340.8(4) 


.02(3) 


-340.8(4) 


-340.6(4) 


q 


-107.6(4) 


.12(4) 


-107.5(4) 


-106.9(4) 


g 
v t 


-4.22(1) - 


.043(1) 


-4.23(1) 


-4.23(1) 


q 


-2.06(1) 


.002(1) 


-2.06(1) 


-2.07(1) 


7T,SR 

Var 


-300(1) 


.18(8) 


-298(1) 


-297(1) 


TT,LR 


23.14(4) 


-.01(1) 


23.13(4) 


23.10(3) 




-9.6(1) - 


.007(2) 


-9.7(1) 


-9.49(2) 



TABLE VIII: Six-quark state energies for cavity diameter R c = 6 fm. Notation identical to Table 
VT 



c s (MeV) A s 


(fm" 1 ) 


2 fm 


4 fm 


6 fm 







8.8(4) 


-1.7(4) 


-1.3(4) 


-0.0077 


5 


-14.2(4) 


-6.5(4) 


-2.6(4) 


-0.67 


2 


-14.2(4) 


-11.0(4) 


-5.1(4) 


&E em p 




-14.2 


-10.5 


-6.0 



TABLE IX: The 5E(R C ) for various c and Ag compared with 5E emp (R c ) from Argonne vis 
potential. 
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0-20 1 




(n) 



FIG. 5: Flux tube configurations for system of six quarks consistent with gauge invariance. Fig. (a) 
shows the flux tubes for the two nucleon state, while (b) shows the "exotic" six-quark configuration. 
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FIG. 6: Flux tube configurations for quark exchange. The solid lines show flux tubes for the 
first partition (123;456). The dashed lines show the flux configuration in the seventh partition 
(163;452), obtained by exchanging quarks 2 and 6. 
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FIG. 7: Internal pair correlation functions in nucleon: (a) central, projected into T, S channels, 
/r 5' (b) S 1 = 1 tensor, // T ; (c) central correlation, (d) operator correlations, = fp/fl, at 
r < 2. They are constant for r > 2 fm. 
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0.5 1 1.5 2 2.5 3 3.5 




0.5 1 1.5 2 2.5 3 3.5 



r(fm) 

FIG. 8: External pair correlation functions in six-quark wave function for R c = 2 fm with 
d c = 3.55, dt = 1.25: central correlation, f^, (top); operator correlations, (bottom). 
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FIG. 11: External central correlation operators projected into T,S channels. 
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FIG. 12: Quark pair distribution functions in six-quark and single nucleon states. The of the 
nucleon is multiplied by a factor of two for easier comparision. 
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I I I I I I I I I I I I I I I I I I I I I 1_ 

0.5 1 1.5 2 

r(fm) 

FIG. 13: The quark pair T q > • T q distribution function in R c = 2, 4, and 6 fm six-quark states and 
in the nucleon. The nucleon p g 2 {r) is multiplied by two for easier comparison. This figure has four 
essentially overlapping curves. 
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